Precision measurements of 5-wave scattering lengths in a two-component 

Bose-Einstein condensate 



M. Egorov/'0 B. Opanchuk,^ P. Drummond,^ B.V. Hall,^ P. Hannaford/ and A. 1. Sidorov^'lll 

^ARC Centre of Excellence for Quantum- Atom Optics and Centre for Atom Optics and Ultrafast Spectroscopy, 
Swinburne University of Technology, Melbourne 3122, Australia 
(Dated: March 2, 2013) 

We use collective oscillations of a two- component Bose-Einstein condensate (2CBEC) of 
^^Rb atoms prepared in the internal states |1) = = mF = — 1) and |2) = |F = 2, mp = ^) 
for the precision measurement of the interspecies scattering length a 12 with a relative uncertainty 
of 1.6 X 10"^. We show that in a cigar-shaped trap the three-dimensional (3D) dynamics of a 
component with a small relative population can be conveniently described by a one-dimensional 
(ID) Schrodinger equation for an effective harmonic oscillator. The frequency of the collective 
oscillations is defined by the axial trap frequency and the ratio ai2/aii, where an is the intra- 
species scattering length of a highly populated component 1, and is largely decoupled from the 
scattering length 022, the total atom number and loss terms. By fitting numerical simulations of 
the coupled Gross-Pitaevskii equations to the recorded temporal evolution of the axial width we 
obtain the value ai2 = 98.006(16) ao, where ao is the Bohr radius. Our reported value is in a 
reasonable agreement with the theoretical prediction a 12 = 98.13(10) ao but deviates significantly 
from the previously measured value ai2 — 97.66 ao [l[ which is commonly used in the characterisa- 
tion of spin dynamics in degenerate ^^Rb atoms. Using Ramsey interferometry of the 2CBEC we 
measure the scattering length a22 = 95.44(7) ao which also deviates from the previously reported 
value a22 — 95.0 ao W\. We characterise two-body losses for the component 2 and obtain the loss 
coefficients 712 = 1.51(18) x 10"^^ cmVs and 722 = 8.1(3) x 10"^^ cmVs. 

PACS numbers: 67.85.De, 67.85.Fg, 34.50.-s, 03.75.Dg 



I. INTRODUCTION 

Collisional interactions in dilute ultracold gases play 
an important role in the dynamics of Bose-Einstein con- 
densates (BECs), formation of molecules and shifts of 
resonance frequencies. Binary collisions of atoms can be 
divided into inelastic and elastic collisions. Inelastic col- 
lisions lead to the change of a hyperfine state or spin flips 
of the colliding atoms. At sufficiently low temperatures 
the elastic collisions of ultracold bosonic atoms can be 
adequately described by a single parameter, the s-wave 
scattering length a, through the corresponding interac- 
tion strength [2]. 

Precise knowledge of the values of the scattering length 
is required for reliable modeling of BEC dynamics and 
spin squeezing, the accurate evaluation of collisional 
shifts in atomic clocks, and for spin gradient thermom- 
etry at sub-nanokelvin temperatures [3]. It can also be 
used to verify theoretical models of inter-atomic poten- 
tials [1]. In general, it is difficult to carry out precision 
measurements of a scattering length in a single species 
ensemble. Binary mixtures in the form of either two con- 
densates of different atomic species or two-component 
Bose-Einstein condensates (2CBEC) provide more op- 
portunities for the accurate measurement of these col- 
lisional properties. Two-component BECs are defined as 
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a mixture of two different spin- or hyperfine states of the 
same condensed species. In either case collisions involve 
intra-species (an and a22) and inter-species (ai2) s-wave 
scattering. 

Due to the particular properties of the singlet and 
triplet inter-atomic potentials [5] the s-wave scattering 
lengths of ^^Rb atoms in the ground hyperfine states with 
F = 1 and F = 2 are very close to each other (the maxi- 
mum difference is around 5%). These two states in ^'^Rb, 
|1) = |F = l,mF = -l) and |2) = |F = 2,mF = l), 
are potentially useful for future applications of ultra- 
cold or condensed atoms because they are magneti- 
cally trappable and their differential first-order Zeeman 
shift is cancelled at a rather low magnetic field strength 
3.228 G [H. As a result, coherent superpositions of these 
two states are largely insensitive to magnetic field noise. 
Under appropriate conditions very long coherence times 
have been reported for trapped non-condensed atomic 
ensembles [7] and a 2CBEC [8]. This makes these states 
of great interest for on-chip atomic clocks and interfero- 
metric applications, since the collisional frequency shift 
in a trapped atomic clock with equal population of the 
two states is proportional to the difference (an — a22) be- 
tween the intra-species scattering lengths. The coherent 
superposition of these states and subsequent nonlinear 
evolution were recently used in a spin-squeezing experi- 
ment [il. 

In this paper we demonstrate a new method for pre- 
cision measurement of the interspecies scattering length 
ai2 using collective oscillations in a 2CBEC. It was pre- 
viously proposed [lo| that the inter-species coupling has 
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a dramatic effect on the collective excitation spectrum. 
Our method is largely decoupled from the a22 scattering 
length, the total atom number and loss terms. In the 
course of our study we measured the two-body loss co- 
efficients for the states |1) and |2) in ^^Rb and obtained 
values that are significantly different from previous mea- 
surements [l|, [ll|. We also carried out a measurement 
of the intra-species scattering length a22 using Ramsey 
interferometry of a trapped 2CBEC, a known theoretical 
value of an and our measured value of ai2. 

In Sec. ini we use the quantum least-action principle 
for a 2CBEC trapped in a cigar-shaped harmonic po- 
tential to derive a ID Schrodinger equation describing 
collective oscillations of the component 2. The oscilla- 
tion frequency is defined by the axial trap frequency and 
the ratio of the scattering lengths ai2/aii. In Sec. HIH 
we use a full three-dimensional simulation of the cou- 
pled Gross-Pitaevskii equations (GPE) with loss terms 
to confirm that the frequency of the collective oscilla- 
tions in the approximation of a small relative popula- 
tion of the component 2 {N2 <^ Ni)^ largely depends on 
the value of ai2 (relative to an) and has a very weak 
dependence on the a22 value and the total atom num- 
ber AT = A^i + 7V2 (for TV > 3 X 10"^). We describe our 
experimental setup and the characterization of the trap 
frequencies in Sec. IIVI Section |Vl describes the converg- 
ing analysis sequence which we use to obtain values of 
the s-wave scattering lengths ai2, a22 and the two-body 
loss coefficients 712 and 722- In Sec. IVIi we present the 
results of our measurement of the two-body loss coef- 
ficients. We describe in detail our measurements of the 
scattering length ai2 using collective oscillations in a two- 
component BEG and a22 using Ramsey interferometry 
in Sections IVIII and IVIII) respectively. We compare our 
results with previous experimental and theoretical inves- 
tigations in Sec. [1X1 



II. ONE-DIMENSIONAL EFFECTIVE 
SINGLE-COMPONENT TREATMENT 

We consider the dynamics of a two-component Bose- 
Einstein condensate initially prepared in the inter- 
nal quantum state |1) (component 1) and trapped in 
a cigar-shaped, axially symmetric harmonic potential 
V = + ma;^r^/2, where = y'^. Electro- 

magnetic radiation can transfer a variable portion of the 
condensate to another internal quantum state |2) (com- 
ponent 2) which is also trapped in the same potential V. 
The transfer modifies mean-field interactions and initi- 
ates a dynamical evolution of the two condensate wave 
functions ^i(r, t) and ^2(1*,^) normalized to the atom 
numbers Ni and N2 in the components 1 and 2 respec- 
tively. In this section, we obtain an analytic expression 
which substantiates decoupling of the ai2 measurement 
from the parameters a22 and N. 

We consider the case of a tight transverse confinement 
(ujr ^ ^z) where the 3D dynamics of the 2GBEG can 



be conveniently described by an effective ID treatment. 
We use the variational method and follow the procedure 
developed for single-component condensates |12l - [l4. The 
action functional of a 2GBEG can be written as [151 
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where the Lagrangian density of the component |j) is 
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where Uj^ 



Aiih^aij/m are the inter- and intra- 



component interaction strengths, and a^j are the s- 
wave scattering lengths (i,j = 1,2). Coupled three- 
dimensional GPE can be obtained as dS/d^^'j = [2!]: 
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In order to reduce the 3D treatment to the ID case, 
we factorize the wavefunctions in the form [l2| 
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where fj is normalized to the atom number in component 
j and (pj is a Gaussian trial function normalized to unity 
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The use of Gaussian trial functions for the radial 
dependence of the condensate density is justified in 
the limit of weak interactions when the BEG is one- 
dimensional [m, [l6| . A typical Thomas- Fermi (TF) ra- 
dius for the BEG in our experiments along the tight trap 
direction is 4 /im, four times larger than the size of the 
corresponding harmonic oscillator ground state. However 
a Gaussian trial wavefunction is known to give consistent 
results for a ID reduction even in the case of a TF ra- 
dial profile of a BEG [12|. We assume that is slowly 



varying along the axial coordinate relative to the radial 
direction and 



Using the Euler-Lagrange equations {dS/dfj 
dS/d(jj = 0) we obtain 
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For the case of the transfer of a smah atom number N2 
to state 1 2) (I/2I <C |/i|) and a large atom number Ni in 
state |1) {\fif > (2aii)-i) we find 
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In the approximation | /2 1 ^ / | /i | ^ <C 1 , the density in 
component |1) does not change and is given in the TF 
approximation by 
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where ji is the effective ID chemical potential of com- 
ponent |1). Equation [13] is valid for |z| < Ttf = 
(2/i/ma;^)^/^ otherwise |/i| = 0. We substitute Eq. [13] 
instead of |/i| for simplicity. However this implies a 
certain limitation on Ttf for which the analytical solu- 
tion is valid, which is discussed after the solution is ob- 
tained (Eq. [21]). Now the effective ID equation for com- 
ponent 1 2) is 
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Apart from the constant term /ieff^ this is the Schrodinger 
equation for a harmonic oscillator. If the superposition 



of states |1) and |2) is prepared by a pulse with area 6>, 
the ID wavefunction of state |2) is expressed in terms of 
a TF profile for state |1) as /2(^,0) = sin^(l9/2) /i(z, 0). 
Therefore, the solution of Eq.[T7] takes the form 



/2(^,t) = e-*^^^^'/^^ 



-iu;eff(2/c+ J) t 



X / ipi.oi'^k.z) f2{z,0)dz 
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where only even harmonic oscillator eigenstates 
?/^ho(2/c,z) contribute to the solution because the wave- 
function 72(^,0) is symmetric about z = 0. Equation [18] 
is periodic in such a way that 



f2{z,t + n/Q = e-^''-»/(^/=)/2(^,t), 



n G 



(19) 




where fc = 2x a;eff/27r gives the frequency of the collec- 
tive oscillations 



(20) 



The effective harmonic potential in equation [15] acts 
only within the size of the BEC whose density is non- 
zero when \z\ < r^F- Therefore, Eq.[T5]is valid only when 
the characteristic size of the relevant harmonic oscillator 
eigenstates is less than Ttf- As a criterion, we require N 
to be large enough to make the first two even eigenstates 
of the harmonic oscillator smaller than Ttf, so that 
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In our experiments the critical value of for this crite- 
rion is 5 X 10^, and we choose A^ to be at least 10 times 
larger than this value. 

Thus the transfer of a small atom number from state 
|1) to state 1 2) initiates collective oscillations of compo- 
nent 2 along the axial coordinate z with frequency fc 
which is independent of the total atom number A" (pro- 
viding criterion[22]holds), the intra-component scattering 
length a22 and the transferred fraction. The evolution 
of the ID density of component 2 (712 = 72/2) is pre- 
sented in Fig. 1(a) and clearly shows periodic compres- 
sions of component 2 with frequency 2.91 Hz (ai2/aii = 
98.006/100.40, = 11.507 Hz). Accurate measurements 
of the axial trap frequency and the frequency of the col- 
lective oscillations make precision measurements of the 
ratio ai2/aii possible. When ai2 < an component 2 
contracts periodically. For ai2 > the dynamics of 
component 2 becomes unstable. 
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FIG. 1. Oscillations of the one- dimension al densit y 1/2(2;, t)|^ 
with frequency 2.91 Hz evaluated from Eg. 1 18) (a7] and the lin- 
ear density with f reque ncy 3.00 Hz simulated with the coupled 
3D GPE (Eqs. ^ [(b)] 



III. GPE SIMULATIONS OF COLLECTIVE 
OSCILLATIONS OF COMPONENT 2 

In this Section we test the accuracy of the predictions 
of the effective ID treatment, by comparing the above re- 
sults with the full numerical simulations of the 3D GPE 
including collisional losses, described by the following 
equations: 
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(23) 



Here the loss rates of species 1 and 2 are Fi = 1712 |^2|^ 

and r2 = f (712 1^1 1^ + 722 1^2!^), and 712 and 722 are 
the two-body loss coefficients [jy|. The three-body loss 
rate is negligible at our typical BEG densities. The 
simulations are performed using a symmetric split-step 
Fourier method [13 • The main idea of this method 



is the integration of the GPE by the separate applica- 
tion of nonlinear and differential operators (the latter 
being applied in Fourier space). The symmetric ver- 
sion of the method additionally applies a differential 
operator in two steps, separated by the application of 
a nonlinear operator in a middle step, allowing it to 
reach a local error of the third order in the time step. 
Free expansion of the BEG is simulated when needed 
on a grid with increased size (128 x 256 x 256 for a 
free expansion compared to 128 x 32 x 32 in a trap) 
with the trapping potential set to 0. We also find the 
optimal regimes for precision measurements of ai2 and 
the dependence on various parameters. Unless other- 
wise specified we use the ^^Rb parameters for the states 
|1) = |F = 1, mF = -1) and |2) = |F = 2, = +1) 
(ail = 100.4 ao, ai2 = 98.006 ao and a22 = 95.44 ao) and 
the trap frequencies (98.23, 101.0, 11.507) Hz. 

For N = 10^ atoms and a transfer of 2.4% of the atoms 
to state 1 2) the evolution of the simulated linear density of 
component |2) {7121 = / ^2^2dxdy^ Fig. |l(b)"]) r esembles 
that from the effective ID treatment (Fig. |l(a')] ) but ex- 
hibits a few additional distinct features. Firstly, the fre- 
quency of the collective oscillations is 3.00 Hz. Secondly, 
additional periodic variations of the linear density with 
frequency 18.3(3) Hz are clearly visible. These represent 
the monopole compression mode, which is also excited 
due to atomic interactions being altered by the transfer 
of atoms to state |2). Indeed, the frequency of these fast 
oscillations is consistent with the estimated value of the 
lowest monopole mode frequency 2fz = 18.2 Hz for 
very elongated traps [HI . 

The frequency of the slow collective oscillati ons is very 
sensitive to variations in the value of ai2 (Fig. 2(a) ) and 
can be approximated in the range of interest by 

/c (62) = (3.00 - 0.63 62 - 0.063 ^1^2) Hz, (24) 

where ^12 = ai2/ao — 98.0 and the preparation pulse 
area is 9 = tt/IO. The simulated value is higher by 3% 
than the value estimated from Eq. [20] which does not 
account for the dynamics of the BEG in the radial direc- 
tion. A one-dimensional treatment accounting for those 
dynamics can be derived [14]; however the resulting equa- 
tions can be solved only numerically. The difference be- 
tween Eq. [20] and the GPE simulations depends weakly 
on N which appears in our simulations as the derivative 
dfc/dN ~ 2 X 10~^. The dependence on a22 is greatly 
suppressed for small atom number in state |2) (Fig. |2(b)"] ). 
For a mixture of two components prepared by a tt/IO- 
pulse, /c can be estimated from 



/c (62) = (3.00 - 0.020 62) Hz, 



(25) 



where £^22 = 0^22/ clq — 95.5. Sensitivity to the total 
atom number TV is also suppressed for small pulse areas 
6 (Fig. 2(c)[ ) and sufficiently large atom numbers (Eq.[22l 
Fig. 2(c)[ ). We have checked that the dependence of 
the frequency on the pulse area is not pronounced when 
6> <C 7r/2. We use the first order derivatives of the collec- 
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FIG. 2. Dependence of the collective oscilla tion s freque ncy 
of component |2) on the scattering lengths ai2 (a) and 022 (b) 
and on the total atom number N 



Dotted lines are the 
analytical predictions of Eq. [2Ql The three-dimensional GPE 
simulations are represented by the solid lines {0 — tt/IO) and 



dashed lines (6' = 7r/5). N = 10^^ for (a) and [(b) 



tive oscillation frequencies on the experimental parame- 
ters in our error analysis later in the paper. The simu- 
lations showed that for sufficiently large N and small 0^ 
the collective oscillations frequency can be used for pre- 
cision measurements of ai2/aii. We have estimated the 
accuracy of measuring ai2 from the trap frequency fz. 
While /c is independent of /r, Eqs. [20] and [23] assume 
that an axial trap frequency measurement with a pre- 
cision Sfz/ fz = 6 X 10~^ (Fig. [3]) leads to an additional 
uncertainty in the ai2 measurement of Sai2 = 0.003 ao- 



IV. EXPERIMENTAL SETUP 

We generate an almost pure condensate of ^ ''Rb atoms 
in state |1) in a cigar-shaped magnetic trap on an atom 
chip [1^ . A two-photon microwave-radiofrequency (MW- 
RF) transition |20| is used for the fast transfer of a 
variable number of atoms from state |1) to state |2). 
The MW radiation (/mw ^ 6.8 GHz) is applied by 
a half-wave dipole antenna located outside of the vac- 
uum chamber and the RF field (/rf ^ 3.2 MHz) is 




t{s) 



FIG. 3. Dipole oscillations of state |1) BEG along the axial 
direction of the magnetic trap. The measured value of the 
axial trap frequency is fz = 11.507(7) Hz. 



coupled using two side wires on the chip [19|. During 
the transfer the MW field is red-detuned from the in- 
termediate state |F = 2,mF = 0) by 1 MHz. The two- 
photon Rabi frequency of the effective two-level system 
is l]i2/(27r) = 500 Hz and the two-photon detuning is 
A/(27r) = 6.8 Hz. 

We use MW spectroscopy of the transition |1) 
\F = 2,mF = 0) to adjust the magnetic field at the trap 
bottom to the value 3.228(5) G at which the first order 
Zeeman shift between states |1) and |2) is cancelled so 
that the atoms in both states experience almost the same 
trapping potential [6] . For accurate knowledge of the to- 
tal atom number N we detect the atoms in both states 
in the same experimental realization. At the end of the 
cycle, after the cloud is released from the trap, we employ 
adiabatic passage with MW radiation [H, [20| to transfer 
the state |1) atoms to the state |3) = |F = 2,mF = — 1) 
with an efficiency of 98%. During the time-of- flight fall 
the atoms in states |2) and |3) are spatially separated 
by a magnetic field gradient generated with a pulse of 
current through the Z-wire on the chip. A single ab- 
sorption image of the atoms in the two states is taken 
using a 100 |Lts pulse of probing light resonant with the 
F = 2 ^ = 3 cyclic transition. The optical reso- 
lution of our imaging system is estimated to be around 
6 |Ltm. The experimental images are post-processed by a 
fringe-removal "eigenface" algorithm [2l| to improve the 
signal-to-noise ratio. 

Accurate knowledge of the harmonic trap frequencies 
(especially fz) is essential for precision measurements 
of the ai2 value. We employ a standard method of 
dipole oscillations by shifting the trap by half of a cy- 
cle, returning to its original position and monitoring the 
periodic oscillations of the BEG in state |1) (Fig. [3]). 
The dipole oscillations in the axial direction are not 
damped over a long period of time. From these measure- 
ments we infer the three harmonic oscillator frequencies 
of our trap: fz = 11.507(7) Hz, = 98.23(5) Hz and 
fy = 101.0(5) Hz. We checked numerically that anhar- 
monicity of the trapping potentials on the atom chip does 
not affect the results. 
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FIG. 4. Iterative convergence of values of the inter-species 
scattering length a 12 



(a) 



and the intra-species scattering 
length a22 |(b)[ Solid lines represent the results obtained 
with a preparation pulse area = t t/IQ , dashed lines are 
for ^ = 7r/5 in (a) Different colors in (a) represent the data 
for different sets of measurements of ai2. 



V. MEASUREMENT SEQUENCE AND 
CONVERGENCE OF ANALYSIS 

The time dependence of the collective oscillations of 
the lightly populated component |2) is crucially depen- 
dent on the ai2/aii ratio. In order to account for the 
residual dependence on the scattering length a22, the to- 
tal atom number N and the two-body loss coefficients, 
we carry out measurements and analysis of the data in 
the following way. We use the theoretical predictions 
of the scattering length values (ai2 = 98.13 ao and a22 = 
95.68 ao) [4] and experimental values of the two-body loss 
coefficients measured at 8.32 G (712 = 7.80 x 10~^^ m^/s 
and 722 = 1.194 x 10~^^ m^/s) [l[ as initial parameters. 
In the first iteration we use these values in the interfero- 
metric calibration of the total atom number N [8[ which 
we find to be consistent with the calibration using the 
condensation temperature. We measure the two-body 
loss coefficients using the atom number calibration re- 
sults. Then we find a new value of ai2 from collective 
oscillations dynamics. In the next step we find a22 from 
Ramsey interferometry measurements with tt/IO and 7r/2 
preparation pulses. At the end of the first iteration we 
find new values of 712, 722, <^12 and a22 and cycle 
through the same sequence of analysis several times until 
all values converge. We find that 3 iterations are suffi- 
cient for convergence (Fig. 2]). 



VI. TWO-BODY LOSS COEFFICIENTS 

In the mean field approximation atom losses of a two- 
component BEG are described by the equations [ll| 



dni 
~dt 
dn2 
~dt 



-712^1^2, 



-722^2 - 712^1^2, 



(26) 



where ni and n2 are the densities of each BEG com- 
ponent, and 7i2 and 722 are two-body loss coeffi- 
cients. During an inelastic collision two condensed 



atoms in state |F = 2,mi? = 1) change their spin states 
to \F = 2,mi? = 0) and \F = 2,mi? = 2). An atom in 
|F = 2, = 0) state is lost from the magnetic trap and 
an atom in \F = 2, mi? = 2) acquires a potential energy 
of 624 nK due to the gravitational sag and moves out of 
the trap. The atoms in these states do not contribute to 
losses of the two-component BEG as they do not over- 
lap with the condensed atoms. Two atoms in states 
\F — 2, mi? = 1) and \F = l,mi? = —1) can spin flip to 
untrappable states |F = 2, ttif = 0) and |F = 1, rriF = 0) 
and do not further contribute to atom losses. Here we 
neglect inelastic collisions with the background gas and 
three-body losses which do not contribute on a timescale 
of less than one second. 

In order to measure the 722 coefficient we produce a 
BEG in state |1) and then prepare a pure state |2) con- 
densate with a TT pulse. After a variable evolution time 

we release the condensate and measure the remain- 
ing number of atoms in state |2). We apply a magnetic 
field gradient to separate in free fall the atoms in |2) and 
\F = 2, mi? = 2) states and only measure the population 
of state 1 2). If the BEG adiabatically follows a TF profile 
during the lossy evolution, the loss of a single component 
is described for short times by [ll| 



^2 Ht) 
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where / = 48.5 Hz is the mean trap frequency. Thus 
the value of 722 is given by the slope of the dependence 
^2"^^^^) (Fig. [5(a)| solid line). In order for the BEG to 
follow the trapping potential adiabatically during the loss 
process, the characteristic loss rate r2 should be less than 
the trapping frequency fz- However, in our experiments 
these quantities are comparable, and we find the atom 
loss by fitting the GPE equations (Eqs. [23]) . However, 
the results of the GPE simulations are very close to the 
linear dependence of Eq. [27] and are not distinguishable 



in Fig. 5(a) 



For the 712 measurements we prepare a superposition 
of two states with a 7r/5 pulse (n2 <C ni) and measure the 
remaining atom numbers in the two components after an 
evolution time t (Fig. |5(b)[ ). Under such conditions the 
loss process depends mostly on 712 rather than 722. This 
allows us to slightly decouple the measurements from the 
value of 722- We find 712 by fitting the experimental data 
with the coupled GPE simulations (Eqs. [23]) using the 
iterated atom number N . 

The exact values of the loss coefficients are strongly 
dependent on the initial total number of atoms (Eq. [27]). 
Using the converging sequence of Section [V] we determine 
the values 712 = 1.51(18) x IQ-^^ m^/s and 722 = 8.1(3) x 
-|^Q-20 j^3^g xhese values differ from the previously re- 
ported values of 722 = 10.4(10) x 10~^^ m^/s measured at 
a magnetic field of 3.0 G in a dipole trap [ll| and 722 = 
11.94(19) X 10-2° mVs and 712 = 7.80(19) x 10-^° m^s 
measured at 8.32 G in a TOP magnetic trap [l|. In 
a separate measurement of the evolution of cold non- 
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FIG. 5. Measurement of intra-species 722 ] (a) [ and inter-species 
712 1(b) I two-body loss coefficients. The plot of N~'^^^ versus 
hold time t is almost linear in the 722 measurement. The 
measurement of 712 is plotted on a semi-log scale, the increase 
in loss rate is observed at the point of maximum density of 
component 2 (0.17 s) due to the influence of 722- Black points 
are experimental results, blue solid lines represent fits with 
GPE simulations. 



condensed atoms prepared with a tt pulse in state |2) we 
determine 722 = 16.9 x 10~^^ m^/s which is by a factor 
of two larger than our measured value of 722- The two- 
particle correlation coefficient for noncondensed bosons 
is equal 2 [g^ , so our measurements of losses of condensed 
and noncondensed atoms are consistent with each other. 



VII. MEASUREMENT OF ai2 SCATTERING 
LENGTH 

We carry out six sets of measurements: three with tt/IO 
pulses and three with 7r/5 pulses. Oscillations of the 
axial width of component 2 are excited (Sec. IIIip and 
we image the column densities of both components in 
time-of- flight after various evolution times (Fig. [6j). The 
experiments are performed with two different times of 
free expansion (6.6 and 20.1 ms) which are also included 
in the GPE simulations. The period of the collective 
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FIG. 6. (a) Absorption images of two BEG components in 
the time-of-flight expansion for different evolution times after 
a 7r/5 preparation pulse, (b-d) Gentral cross sections of the 
column density of component 2 at evolution times of (b), 
100 (c) and 170 (d) ms. Black dots represent the measured 
optical density from the GGD pixels. Blue solid lines are the 
central cuts of the 2D Gaussian functions with axial widths 
of 21.7 (b), 11.1 (c) and 7.6 (d) ^im. 



oscillations depends on the ratio ai2/aii rather than just 
the value of ai2 fE g.lSOl) and we use the established value 
an = 100.40 ao U^U^ 

The two dimensional distribution of the column den- 
sity of the second component is fitted with a 2D Gaus- 
sian function and from here we extract axial widths of 
the column density profiles. The axial cross sections in 
the centre of the 2D profiles are shown in Figs. [6fb-d). 
The choice of a Gaussian function to fit the experimental 
data originates from the fact that Eq. [15] assumes that 
the ground state of component 2 has a Gaussian shape. 
Even when the initial BEG density profile has the shape 
given by the TF approximation, a Gaussian function fits 
the experimental cross section well (Fig. [Gb) and is a 
good measure of the BEG width. In our analysis we do 
not use second moments of the column density for width 
measurements because in this case the extracted BEG 
width has error bars larger by a factor of ~ 20. 

Oscillations of the axial width of component 2 with 
time are shown on Figs. [7(a)] and [7(b)] for tt/IO and 7r/5 
preparation pulses and a time-of-flight of 6.6 ms. Each 
point represents one experimental realization and error 
bars are the statistical uncertainty of the Gaussian fits 
of the recorded 2D column densities. We fit the tempo- 
ral dependence of the axial width with the coupled GPE 
equations (Eq. [23|) varying ai2, a22 and N and using the 
iteration procedure of Sec. [Vj The extracted values of 
the scattering length ai2 for the six measurement sets are 
shown in the Table [T] Systematic errors are calculated for 
each set of measurements assuming 10% uncertainty in 
the preparation pulse area 6 and 15% preparation noise 
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FIG. 7. Temporal evolution of the axial width of com pone nt 
2 in a superposition of two states prepared with a 7r/lQ |(a)] or 
7r/5 |(b)| pulse. The expansion time is 6.6 ms and the initial 
total atom number N = 1.1 x 10^. Black dots are the data 
extracted with the 2D Gaussian functions and error bars rep- 
resent statistical errors of the fits. Blue solid lines are the re- 
sults of the GPE simulations (Eq.[23l) with ai2 = 97.961 ao (a) 
and 97.986 ao (b), and 022 = 95.44 ao. 
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FIG. 8. Ramsey interferometry with a variable area for the 
preparation pulse is used in measurements of the scattering 
length a22- |(a)| The normalized number difference Pz oscil- 
lates with the evolution time t. Dots - data points, blue solid 
line - fitted GPE simulations for — tt/IO, red solid line - fit- 
ted GPE simulations for = tt/2. |(b)| Two-photon detuning 
A/(27r) — 6.82(10) Hz was measured in Ramsey interferome- 
try of non-condensed atoms. Dots - data points, solid line - 
sinusoidal fit. 



dependencies (Fig. [2j) 
Sa 

\(Jc 
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The final value of ai2 = (98.006 ± 0.016) ao is calcu- 
lated as the weighted mean of the six results, where the 
weight coefficients are obtained from the errors of indi- 
vidual measurements (Table H]). The inclusion of quantum 
noise [23[ and finite temperature effects [2^ can improve 
the precision of the measurements and will be the subject 
of further studies. 
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Weighted mean: ai2 (98.006 ± 0.016) ao 



TABLE I. Values of the scattering length ai2 extracted from 
the six experimental sets with uncertainties coming from dif- 
ferent sources. 



in the initial total atom number N. The combined un- 
certainty of ai2 consists of the fit error Siai2 and the sys- 
tematics calculated from the slopes of the corresponding 



VIII. MEASUREMENT OF a22 SCATTERING 
LENGTH 

In this series of measurements we perform Ramsey 
interferometry of a BEG initially prepared in state |1) 
using di = 7r/2 (or = tt/IO) preparation pulse 
and a 7r/2 interrogation pulse (Fig. [8] top). Our mea- 
sured quantity, the normalized number difference Pz = 
{N2-Ni)/{N2^Ni) (Figs. [8(a3 and [(b)}, where TVi and 
N2 are the populations of the two states after the interro- 
gation pulse, oscillates with frequency determined by the 
detuning of the MW/RF radiation from the two-photon 
resonance and the collisional shift. Two experimental se- 
quences are intermixed in time to ensure that the data 
points are taken for the same total atom number N and 
two-photon detuning. For the detuning measurement we 
perform Ramsey interferometry (two 7r/2 pulses) with a 
non-condensed cold atomic ensemble {N ^ 1.8 x 10^, 
T ^ 120 nK, Fig. |8(b)[ ) in the same magnetic trap and 
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find the two-photon detuning A = 6.82(10) Hz. For 
the cold atoms the collisional shift is estimated to be 
- 0.1 Hz. 

The collisional shift of dense condensed atoms de- 
pends not only on the particular values of the scatter- 
ing lengths ail, <^12 and a22 but also on the total num- 
ber of atoms TV. We use two different values of the 
preparation pulse area for the following reason. If 
we assume that the atom number density n(r) does not 
change with time the collisional shift is proportional to 
n(r)(aii — a22) for 6> = 7r/2 [6[. For a 6> = tt/IO prepa- 
ration pulse the interaction-induced shift is proportional 
to n(r)(aii — ai2). Thus the ratio of the two shifts is 
independent of the atom number density and, therefore, 
independent of TV. In order to obtain a precise value of 
a22, we run the GPE simulations with iterative values 
of the scattering lengths (Sec t. |Vj) and fit the Ramsey 
fringe of = tt/IO (Fig. 8(a)[ solid line), keeping N as 
a free parameter and taking into account that the fringe 
frequency is largely decoupled from the a22 value. Us- 
ing the fitted value of N we fit the ^ = 7r/2 data with 
the GPE simulations with the single variable parameter 
a22 (Fig. 8(a)[ dashed line). After convergence of the 
analysis of all measurements the value of the intraspecies 
scattering length for atoms in state |2) is found to be 
a22 = 95.44(7) a^. 



IX. DISCUSSION 

A number of different approaches have been used to 
evaluate the s-wave scattering lengths of ^^Rb atoms pre- 
pared in the two lowest hyperfine states. Here we restrict 
our discussion to the modeling and measurements of the 
scattering lengths an, a22 and ai2 at reported values of 
the magnetic field. 

An early calculation of the value an = 106(6) ao (25| 
was based on knowledge of the 5-wave bound-state ener- 
gies of ^''Rb2. Later, using the results of high-resolution 
molecular spectroscopy, the parameters of the atomic 
scattering potentials were calculated which allowed eval- 
uation of an = 100.40(10) ao and ai2 = 98.175 ao [5] at 
a magnetic field of 3.2 G. Cornell's group used the same 
parameters of the scattering potential to evaluj^te an = 
100.44 ao, ai2 98.09 ao and a22 95.47 ao [6] for the 
same magnetic field. Improved modeling of the atomic 
interaction potentials [^] produced the latest theoretical 
values of the scattering lengths: an = 100.40(10) ao, 
ai2 = 98.13(10) ao and a22 = 95.68(10) ao [26]. 

Experiments with two-component BECs provided an- 
other way to measure the scattering properties of 
^^Rb atoms. The transfer of the entire population from 
state |1) to state |2) produced a sudden change of col- 
lisional interactions in the condensate and triggered ra- 
dial and axial collective oscillations of component 2 in a 
time-orbiting-potential trap [l^- GPE simulations were 
used to fit the temporal evolution of the radial and ax- 
ial widths and yielded the ratio aii/a22 = 1.062(12). If 





ai2/ao a22/ao 


Matthews et al [27] 
van Kemp en et al. 
Harber et al. ^ 
Mertes et al. \l\ 
Kokkelmans [26] 


94.5(1.1) 

98.175 

98.09 95.47 
97.66 95.0 
98.13(10) 95.68(10) 


This work 


98.006(16) 95.44(7) 



TABLE II. Calculated and measured values of s-wave scatter- 
ing lengths in ^^Rb for the magnetic field of 3.2 G. Paper [l|] 
reported the results for the magnetic field of 8.32 G. 



we assume an = 100.40 ao then this measurement yields 
a22 = 94.5(1.1) ao. Collisional shift measurements [H 
in Ramsey interferometry of a two-component BEG and 
uncondensed atoms at the bias magnetic field of 3.23 G 
produced a value for the difference of scattering lengths 
<^ii — GL22 = 4.85(31) ao. Collisionally driven spin oscilla- 
tions of atom pairs trapped in an optical lattice at 0.24 G 
allowed measurements of the scattering length differences 
in the F = 1 and 2 manifolds [22|. The reported results 
consistently deviated by up to 10% from the calculated 
values [5]. Use of the computed value an = 100.40 ao and 
a fit of the observed oscillating ring-like structures in a 
two-component BEC at a bias magnetic field of 8.32 G 
with the simulations of the coupled GPE equations have 
yielded ai2 = 97.66 ao and a22 = 95.0 ao [1]. These three 
values of the scattering lengths are commonly used for 
modeling the Ramsey contrast and spin squeezing evolu- 
tion in atom c hip experiments on entangled atomic en- 
sembles carried out at a bias magnetic field of 
3.23 G or spin dynamics of 2CBEC in \F = l,mF = 1) 
and \F = 2,mF = -1) states [13, [sH. 

In Table nil we compare the results of our measurements 
with previously reported values of ai2 and a22. Our value 
of ai2 is within 0.17% of the theoretical values [El. : 
however the statistical uncertainty of our measured value 
is ten times smaller. All reported theoretical values are 
consistently shifted towards higher values. Our value of 
ai2 = 98.006 ao deviates from that of [l] by 0.35%. The 
dependence of ai2 on the magnetic field can explain this 
very significant deviation. Our measured value of a22 
is very close to that of [3] and is within 0.25% of the 
latest theoretical evaluation [26| but is well outside the 
statistical uncertainty 0.07% of our measurement. 

Precise knowledge of collisional loss rates is important 
as this limits the BEC coherence time [1, [2^ , the effec- 
tiveness of spin-squeezing [23,[32[ and contains important 
information about the properties of inter-atomic inter- 
action potentials [sj]. Two-body losses are usually the 
dominant source of losses in binary BECs. 

In this work we have focused on collisions of 
^'^Rb atoms in states |1) and |2) and the two-body loss co- 
efficients 7i2 and 722- The coefficient 712 was initially cal- 
culated work by analyzing spectroscopic studies [5[ . The 
reported value of the imaginary part of a complex scat- 
tering length lm{a{B)) = — 0.02ao at the magnetic field 
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7i2(m^/s) 722(m^/s) 

van Kempen et al [5] 1.9 x 10~^° 

Mertes et al [Ij 7.80(19) x 10"^° 11.94(19) x 10"^° 

Tojo et al [llj 10.4(10) x 10"^° 

This work 1.58(18) x 10"^" 8.1(3) x 10~^° 

TABLE III. Calculated and measured values of two-body loss 
coefficients in ^^Rb for magnetic fields of 3.2 G H, 8.32 G d 
and 3.0 G 11]. 

of 3.23 G can be used to evaluate the two-body loss rate 
by 712 = -4/i/m lm{a{B)) = 1.9 x 10"^^ m^/s. Also this 
work predicted a weak Feshbach resonance at a magnetic 
field of 1.9 G (where 712 increases by factor of 2) which, 
however, has not yet been observed. The two-body loss 
coefficients 712 = 7.80 x 10"^^ m^/s and 722 = 11.94 x 
-|^Q-20 ]^3^g y^QYe first experimentally characterized at the 
magnetic field of 8.32 G [1]. Another experimental mea- 
surement provided 722 = 10.4(10) x 10~^^ m^/s at a bias 
magnetic field of 3.0 G assuming the two-body loss co- 
efficient for state \F = 2,mi? = — 1) to be the same as 
for \F = 2,mF = +1) at low fields jll|. The latter work 
carried out a comprehensive study of two-body losses of 
F = 2 ^^Rb atoms trapped in an optical trap. All mag- 
netic states were trapped and the presence of = or 
+2 atoms may have influenced the reported results. 

In Table IIIII we compare the results of our measure- 
ments with previous theoretical and experimental inves- 
tigations. Our value of the inter-species two-body loss 
coefficient 712 is much closer to the theoretical result [5| 
than the previous experimental result [1]. Also our mea- 
surement of 722 is slightly smaller than the results ob- 
tained in other experimental works [l|, |Tl| . 

X. CONCLUSION 

We have presented a new technique for precision 
measurement of the inter-species scattering length ai2 
in two-component Bose-Einstein condensates which em- 
ploys collective oscillations of a less populated compo- 
nent. The oscillations can be described either by sim- 



ulations using coupled Gross-Pitaevskii equations or by 
the analytical model which we have developed. System- 
atic errors such as uncertainty in the total number of 
atoms, the effect of one of the inter-species scattering 
lengths and imperfect preparation of the two-component 
mixture contribute very little to the uncertainty of the 
measured ai2. We have applied the technique to measure 
the scattering length ai2 = 98.006(16) ao for states |1) = 
\F = 1, mF = -1) and |2) = |F = 2, = 1) in ^'^Rb us- 
ing a calculated value of an = 100.40 which was as- 
sumed to be fixed and error- free. We also have measured 
the intra-species scattering length a22 = 95.44(7) ao us- 
ing Ramsey interferometry and the two-body loss coeffi- 
cients 712 = 1.58(18) X 10-2° m^/s and 722 = 8.1(3) x 
-|^Q-20 j^3y/g fitting atom losses in the two-component 
BEG with the results of simulations of the coupled GPE. 

Comparison of our results for ai2 and a22 with the the- 
oretical predictions and experimental measurements are 
presented in Tab. [Ill Our results show good agreement 
with recent theoretical calculations [26| . The residual de- 
viations from the theoretical predictions could be caused 
either by uncertainty in the theoretical method or by 
unaccounted effects in our experiments. The magnetic 
trapping potential on the atom chip is slightly anhar- 
monic; however we have found that anharmonicity does 
not affect the GPE simulations. Given the high precision 
of our measurements, the results could be affected by 
quantum dynamics beyond mean-field theory. Another 
unaccounted contribution is from the finite temperature 
of the ensemble. We will include quantum dynamics at 
finite temperatures into the simulations in future exper- 
iments [23]. 

The present technique of ai2 measurement might be 
extremely useful in characterization of narrow inter- 
species Feshbach resonances, such as RF-induced Fesh- 
bach resonances [33- 35] . Precision measurements of scat- 
tering lengths in the vicinity of narrow Feshbach reso- 
nances might also allow drifts of electron/proton mass 
ratio to be monitored [s^. 

We acknowledge helpful discussions with S. Whitlock. 
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cellence for Quantum- Atom Optics, and an ARC LIEF 
grant LE0668398. 
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